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Abstract 

We give an overview of the derivation of multipolar equations of 
motion of extended test bodies for a wide set of gravitational theo¬ 
ries beyond the standard general relativistic framework. The classes 
of theories covered range from simple generalizations of General Rel¬ 
ativity, e.g. encompassing additional scalar fields, to theories with ad¬ 
ditional geometrical structures which are needed for the description of 
microstructured matter. Our unified framework even allows to han¬ 
dle theories with nonminimal coupling to matter, and thereby for a 
systematic test of a very broad range of gravitational theories. 


1 Introduction 

In this work we present a general unified multipolar framework, which en¬ 
ables us to derive equations of motion of extended test bodies for a wide 
range of gravitational theories. The framework presented here can be ap¬ 
plied to theories which significantly go beyond General Relativity (GR), and 
range from the most straightforward extensions of GR, like scalar-tensor 
theories, to theories with additional geometrical structures and nonminimal 
coupling. 

The multipolar method which we employ here can be thought of as the 
direct generalization of the ideas pioneered by Mathisson [T], Papapetrou 
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Figure 1: General idea behind multipolar approximation schemes: The 
world-tube S of a body is replaced by a representative world-line L, whereas 
the original energy-momentum tensor is substituted by a set of multipole 
moments ' along this world-line. Such a multipolar description simpli¬ 
fies the equations of motion. This is achieved by consideration of only a finite 
set of moments. Different flavors of multipolar approximation schemes exist 
in the literature, in this work we define the moments a la Dixon in [3]. 


[2], and Dixon [31 HJ [H E] to the case of generalized gravity theories. As 
sketched in figure [H the main aim of such methods is to find a simplified 
local description of the motion of extended test bodies in terms of a suitable 
set of multipolar moments, which catches the essential properties of the 
body at the chosen order of approximation. 

In this work, we use a covariant multipolar description, based on Synge’s 
world-function formalism mM- The multipolar moments to be introduced, 
can be viewed as a direct generalization of the moments first introduced by 
Dixon in |3]. 

Central to the derivation of the equations of motion, by means of a 
multipolar method, is the knowledge of the corresponding conservation laws 
of the underlying gravity theory. In General Relativity, the starting point of 
all methods mentioned so-far is the conservation of the (symmetric) energy- 
momentum 

= 0 . ( 1 ) 

Many generalized gravity theories allow for a more detailed description of 
matter than standard GR, e.g. also taking into account internal matter 
degrees of freedom as source of the gravitational field. Consequently the 
conservation laws of such theories are more complex, and © has to be 
replaced by a suitable generalization. 

Generalized gravity A metric and connection are the two fundamental 
geometrical objects on a spacetime manifold. They play an important role in 
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the description of gravitational phenomena in the framework of what can be 
quite generally called an Einsteinian approach to gravity. The principles of 
equivalence and general coordinate covariance are the cornerstones of this 
approach. As Einstein himself formulated, the crucial achievement of his 
theory was the elimination of the notion of inertial systems as preferred 
ones among all possible coordinate systems. 

In Einstein’s theory, gravitation is associated with the metric tensor 
alone. Nevertheless, it is worthwhile to stress that Einstein clearly un¬ 
derstood the different physical statuses of the metric and the connection 
(Appendix II in [9]): 

“[...] at first Riemannian metric was considered the funda¬ 
mental concept on which the general theory of relativity and thus 
the avoidance of the inertial system were based. Later, however, 
Levi-Civita rightly pointed out that the element of the theory 
that makes it possible to avoid the inertial system is rather the 
inhnitesimal displacement field E . The metric or the symmet¬ 
ric tensor field gik which dehnes it is only indirectly connected 
with the avoidance of the inertial system insofar as it determines 
a displacement field.”. 

There exists a variety of gravitational theories that generalize or extend 
the physical and mathematical structure of GR. Among these theories there 
are large classes of so-called f{R) models, and of theories with nonmini- 
mal coupling to matter; they are developed in particular in the context of 
relativistic cosmology (but not only there), see [TOl dH [El IEl E] ■ The so- 
called Palatini approach represents another class of widely discussed theories 
in which the metric and the connection are treated as independent variables 
in the action principle daiiiiiE]. Another early example of a generalized 
gravity theory with balance laws different from Einstein’s gravity can be 
found in [EllIS]. Last but not least, we should mention the vast family of 
the gauge gravity theories constructed using a Yang-Mills type of approach 

[sniiii]. 

Since our main aim is to present a unified multipolar framework for a 
very wide range of gravity theories, the theory underlying our analysis has to 
be sufficiently general. In this work we choose metric-affine gravity (MAG) 
as the foundation for the derivation of the equations of motion of extended 
deformable test bodies. 

In MAG, matter is characterized by three fundamental Noether currents 
- the canonical energy-momentum current, the canonical hypermomentum 
current, and the metrical energy-momentum current. These objects satisfy a 
set of conservation laws (or, more exactly, balance equations). In view of the 
multi-current characterization of matter in metric-affine gravity, we develop 
a general approach which is applicable to an arbitrary set of conservation 
laws for any number of currents. The latter can include gravitational, elec¬ 
tromagnetic, and other physical currents if they are relevant to the model 
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Figure 2: In contrast to GR, MAG also allows to take into account the 
microstructural properties of matter (spin, dilation current, proper hyper¬ 
charge). These additional currents couple to the post-Riemannian geometric 
degrees of freedom of the underlying spacetime. Our unified description of 
extended test bodies allows to cover a wide range of gravitational theories 
and matter models. 
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under consideration. In particular, our approach is flexible enough to be 
applied to the case in which there is a general nonminimal coupling between 
gravity to matter. Our presentation here is mainly based on the original 
works [22l|23l[Sj. 

Structure of the paper The structure of the paper is as follows: In 
section [2] we give an overview of the metric-affine theory of gravity. In par¬ 
ticular we cover its geometrical and dynamical aspects. This is followed by 
a general Lagrange-Noether analysis in section [3l The results of this gen¬ 
eral analysis are then used to derive the conservation laws of a metric-affine 
gravity with nonminimal coupling in section HI In section [5] a unified multi¬ 
polar framework is presented. This framework is then used [6] to derive the 
general equations of motion of extended test bodies in nonminimal metric- 
affine gravity. Several specializations to other gravity theories are discussed. 
Furthermore, in section [3 the equations of motion in scalar-tensor theories 
are worked out. Section [8] contains a summary of our results. 

Our conventions are those of [25]. In particular, the basic geometrical 
quantities such as the curvature, torsion, and nonmetricity are defined as 
in |25|, and we use the Latin alphabet to label the spacetime coordinate 
indices. Furthermore, the metric has the signature It should 

be noted that our definition of the metrical energy-momentum tensor differs 
by a sign from the definition used in [m 1111126]. A summary of our notation 
and conventions, as well as additional material for comparison to previous 
works can be found in the appendix. 

2 Metric-afRne gravity primer 

Metric-affine gravity |25j is a natural extension of Einstein’s general rela¬ 
tivity theory. It is based on gauge-theoretic principles [MIEI], and it takes 
into account microstructural properties of matter (spin, dilation current, 
proper hypercharge) as possible physical sources of the gravitational field 
on an equal footing with macroscopic properties (energy and momentum) 
of matter. The formalism of MAG makes it possible to study all the above- 
mentioned alternative theories in a unified framework. The corresponding 
spacetime landscape m includes as special cases the geometries of Riemann, 
Riemann-Cartan, Weyl, Weitzenbock, etc. (c.f. figure [3]). 

The standard understanding of metric-affine gravity is based on the 
gauge-theoretic approach for the general affine symmetry group. In this 
framework, we can naturally distinguish the kinematics of the gravity the¬ 
ory and its dynamics. The kinematics embraces all aspects that are related 
to the description of the fundamental variables, their mathematical prop¬ 
erties and physical interpretation. The dynamics studies the choice of the 
Lagrangian and the field equations. We collect in this section all the relevant 
material. Although the Lagrange-Noether machinery is deeply interrelated 
with both kinematics and dynamics of the theory, we will discuss it in a sep- 
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arate section which will ultimately underlie the derivation of the multipole 
equations of motion. 


2.1 Kinematics of metric-afRne gravity 

We model spacetime as a four-dimensional smooth manifold. In this study 
we are not interested in the global (topological) aspects, and we will confine 
our attention only to local issues. The local coordinates x®, i = 0,1,2,3, 
are introduced in the neighborhood of an arbitrary point of the spacetime 
manifold. The geometrical (gravitational) and physical (material) variables 
are then fields of different nature (both tensors and nontensors) over the 
spacetime. They are characterized by their components and transformation 
properties under local diffeomorphisms x® ^ x® + fe®, where 

<5x® = ^®(x). (2) 

The four arbitrary functions ^®(x) parametrize an arbitrary local diffeomor- 
phism. 

2.1.1 Geometrical variables 

In the metric-affine theory of gravity, the gravitational physics is encoded in 
two fields: the metric tensor gij and an independent linear connection 
The latter is not necessarily symmetric and compatible with the metric. 
From the geometrical point of view, the metric introduces lengths and angles 
of vectors, and thereby determines the distances (intervals) between points 
on the spacetime manifold. The connection introduces the notion of parallel 
transport and defines the covariant differentiation Vfc of tensor fields. 

Under the spacetime diffeomorphisms m, these geometrical variables 


transform as 

^Sij = - 9kj - gik-, (3) 

= -{dke)Tu^-{die)Tkf+ {di^i)Tj-dUK (4) 

In general, the geometry of a metric-affine manifold is exhaustively char¬ 
acterized by three tensors: the curvature, the torsion and the nonmetricity. 
They are defined m as follows: 

RkU^ ■= dkVu^ -diTki^ +Tkn^Tu^ (5) 

Tki^ := Vki^-Tik\ ( 6 ) 

Qkij • ^k9ij dk9ij T gij -|- T^j gn. (7) 


The curvature and the torsion tensors determine the commutator of the 
covariant derivatives. For a tensor of arbitrary rank and index 

structure: 


(V.Vfe - = -Tab^^eA^^-^^d^..A 

k I 

i=l j=l 
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Figure 3: Different spacetime types as special cases of a general metric-affine 
spacetime {Rku^ / 0, TkC / 0, Qkij 7 ^ 0). The abbreviations R (curvature), 
T (torsion), and Q (nonmetricity) over the arrows denote the vanishing of 
the corresponding geometrical object. 


The Ricci tensor is introduced by Rij := Rkij^, and the curvature scalar is 
R := Rij. 

A general metric-affine spacetime {Rkit 7 ^ 0, TkC 7 ^ 0, Qkij 7 ^ 0) incor¬ 
porates several other spacetimes as special cases, see figure[3]for an overview. 
The Riemannian connection T^j* is uniquely determined by the conditions 
of vanishing torsion and nonmetricity which yield explicitly 

bfcj ~g {pjQki T ^k9ij di9kj)- (9) 

Here and in the following, a tilde over a symbol denotes a Riemannian 
object (such as the curvature tensor) or a Riemannian operator (such as 
the covariant derivative) constructed from the Christoffel symbols ([9]). The 
deviation of the geometry from the Riemannian one is then conveniently 
described by the distortion tensor 

iVfe/ := ffc/ - Tfc/. (10) 

The system Q and ([7D allows to find the distortion tensor in terms of the 
torsion and nonmetricity. Explicitly: 

Nkj^ = -\{Tkj^ + T\j + Tjk) + - Qkj' - Qjk). (11) 

Conversely, one can use this to express the torsion and nonmetricity tensors 
in terms of the distortion. 


Tfc/ = -2N[kj]\ 

Qkij ‘2Nk(ijy 
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( 12 ) 

(13) 












Substituting m into ©, we find the relation between the non-Riemannian 


and the Riemannian curvature tensors 

Radc’’ = Radc - VaNdJ^ + V dNac^ + Nan^Ndc^ - Ndn^Nac^ . (14) 

Applying the covariant derivative to (HI)-© and antisymmetrizing, we 
derive the Bianchi identities m- 

V[n-RfcZ]/ = T[kl'^Rn]mi, (15) 

(16) 

^[nQk]ij Rnk{ij)‘ (H') 


2.1.2 Tensors, densities, and covariant differential operators 

Along with tensors, an important role in physics is played by densities. A 
fundamental density y/—g is constructed from the determinant of the metric, 
g =deigij. Under diffeomorphisms ([2]) it transforms as 

= - {diO ( 18 ) 

This is a direct consequence of ([3|). From any tensor one can construct 

a density = y/—gBi,,,^'". Although in this paper we will encounter 

only such objects, it is worthwhile to notice that not all densities are of this 
type - see the exhaustive presentation in the book of Synge and Schild [25]. 

There are two kinds of covariant differential operators on the spacetime 
manifold, depending on whether the connection is involved or not. The 
Lie derivative is defined along any vector held C* it maps tensors 
(densities) into tensors (densities) of the same rank. Let us recall the explicit 
form of the Lie derivative of the metric and the distortion: 


Rc9ij = C^dkgij + {diC^)gkj + {djC,^)gik, (19) 

= CdnNkj^ + {dkCWnj^ + {djC)Nkn'-{dnC)Nkj^. (20) 

In contrast, a covariant derivative raises the rank of tensors (densi¬ 
ties) and it is determined by the linear connection F^j*. Moreover, there are 
different covariant derivatives which arise for different connections that may 
coexist on the same manifold. 

A mathematical fact is helpful in this respect: Every third rank tensor 
Afcj* dehnes a map of one connection into a different new connection 

Ffc/ ^ Ffc/ + Afc/. (21) 

There are important special cases of such a map. One example is obtained 
for Xkj^ = Nkj^'- then the connection F^j* is mapped into the Riemannian 
Christoffel symbols, F^/ = F^j* -|- Nkj"^, in accordance with ([TO]) . 

Another interesting case arises for Xkj^ = Tjk. The result of the map- 
ping _ 

^ kj — kj “r -L jk — ^ kj ' ^ jk kj — ^ jk • 
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is then called a transposed connection, or associated connection, see [2^ 150] . 

The importance of the transposed connection is manifest in the following 
observation. Although the Lie derivative is a covariant operator - this is not 
apparent since it is based on partial derivatives - one can make everything 
explicitly covariant by noticing that it is possible to recast and (120]) 
into equivalent form 

kdij + {^iC^)gkj + iy jC^)9ik, (23) 

yNki = CVnNk/ + {VkONni + yjC)NkJ - (24) 

By the same token we can “covariantize” the Lie derivatives for all other 
tensors of any structure and of arbitrary rank. 

A more nontrivial (and less known) fact is that we can define the Lie 
derivatives also for objects which are not tensors. In particular, the Lie 
derivative of the connection then reads [29] : 

yTki = ^kVjC - Rki/C'- (25) 

This quantity measures the noncommutativity of the Lie derivative with the 
covariant derivative 


k I 

i=l j-1 

The connection Tj./, the transposed connection and the Rieman- 
nian connection Vy define the three respective covariant derivatives: V^, 
V k, and V k ■ The covariant derivative defined by the Riemannian connection 
is conventionally denoted by a semicolon . 

We will assume that these differential operators act on tensors. In ad¬ 
dition, we will need the covariant operators that act on densities. For an 
arbitrary tensor density ' we introduce the covariant derivative 


:= (27) 

which produces a tensor density of the same weight. We denote a similar dif¬ 
ferential operation constructed with the help of the Riemannian connection 

by 

:= (28) 

When we find 

(29) 

(30) 

where we introduced a modified covariant derivative 

Vi-=Vi + Nky (31) 
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2.1.3 Matter variables 

We will not specialize the discussion of matter to any particular physical 
field. It will be more convenient to describe matter by a generalized field 
The range of the indices A,B,... is not important in our study. How¬ 
ever, we do need to know the behavior of the matter field under spacetime 
diffeomorphisms ([2]): 

= (32) 

Here (ct^b)/ are the generators of general coordinate transformations that 
satisfy the commutation relations 

= {cj^bV 5^ - 5?. (33) 

We immediately recognize in (I33p the Lie algebra of the general linear group 
GL(4, ii). This fact is closely related to the standard gauge-theoretic in¬ 
terpretation [25] of metric-affine gravity as the gauge theory of the general 
affine group GH(4, i?), which is a semidirect product of spacetime translation 
group times GL(4,12). 

The transformation properties (I32p determine the form of the covariant 
and the Lie derivative of a matter field: 


:= - Tfci^ (a^s)/V’^, (34) 

:= + (5iC^)(<7^s)i* (35) 

The commutators of these differential operators read 

(VfcV; - V;Vfc)V’^ = -RMj\o^B)iH^ (36) 

(T^Va: - = -{C(Tki){a^B)iH^. (37) 


2.1.4 Symmetries in MAG: generalized Killing vectors 

As is well known, symmetries of a Riemannian spacetime are generated 
by Killing vector fields. Each such field defines a so-called motion of the 
spacetime manifold, that is a diffeomorphism which preserves the metric 
9ij ■ 

Suppose C is a Killing vector field. By definition, it satisfies 

y iCj ~^y jCi — 0 - ( 38 ) 

By differentiation, we derive from this the second covariant derivative 

VjVjCfe = RjkiCi- (39) 

Apply another covariant derivative and antisymmetrize: 

V[„Vj]VjCfc = y[n{R\jk\i]’'Ci)- (40) 
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After some algebra, the last equation is recast into 

riRijkl + Rnjkl^iC^ + R-inkl^ jC^ + Rijnl^ kC^ + RijkvS iC^ = 0- (41) 
The equations ([381) . (p9]l and (|1T]) have a geometrical meaning: 

R(^9ij — 0, (42) 

= 0, (43) 

R-C^ijki = 0. (44) 

That is, the Lie derivatives along the Killing vector field C vanish for all 
Riemannian geometrical objects. Moreover, one can show that the same is 
true for all higher covariant derivatives of the Riemannian curvature tensor 

m 

R-C . . . Vrif^Rijkl^ = 0- (45) 

Let us generalize the notion of a symmetry to the metric-affine spacetime. 
We take an ordinary Killing vector field ^ and postulate the vanishing of the 
Lie derivative 

= 0 (46) 

of the distortion tensor. Combining this with dill) and d43l) . we find an 
equivalent formulation 


R-Cdij — 0) (47) 

= 0. (48) 

We call a vector field that satisfies (j47)) and (f48]) a generalized Killing vector 
of the metric-affine spacetime. By definition, such a C, generates a diffeo- 
morphism of the spacetime manifold that is simultaneously an isometry (j47[) 
and an isoparallelism d48l) . 

Since the Lie derivative along a Killing vector commutes with the covari¬ 
ant derivative, see (1261) . we conclude from (fT4]) and (1441) that 

the generalized Killing vector leaves the non-Riemannian curvature tensor 
invariant 

C^Rkij' = 0. (49) 

It is also straightforward to verify that 

CC (Vni . . . ^n^Rkli) = 0 (50) 

for any number of covariant derivatives of the curvature. 

Later we will show that generalized Killing vectors have an important 
property: they induce conserved quantities on the metric-affine spacetime. 
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2.2 Dynamics of metric-afRne gravity 

The explicit form of the dynamical equations of the gravitational field is 
irrelevant for the conservation laws that will form the basis for the deriva¬ 
tion of the test body equations of motion. However, for completeness, we 
discuss here the field equations of a general metric-affine theory of gravity. 
The standard understanding of MAG is its interpretation as a gauge theory 
based on the general affine group GA(4, ii), which is a semidirect product 
of the general linear group GL(4, i?), and the group of local translations 
|25j . The corresponding gauge-theoretic formalism generalizes the approach 
of Sciama and Kibble [32l|33]; for more details about gauge gravity theories, 
see [2D1 El] • In the standard formulation of MAG as a gauge theory [2S] , 
the gravitational gauge potentials are identified with the metric, coframe, 
and the linear connection. The corresponding gravitational field strengths 
are then the nonmetricity, the torsion, and the curvature, respectively. 

Here we use an alternative formulation of MAG, in which gravity is de¬ 
scribed by a different set of fundamental field variables, i.e. the independent 
metric gij and connection T ki ^, see [Ml EH EH EH EH E3 EH EH EH [40] , 
for example. It is worthwhile to compare the field equations in the differ¬ 
ent formulations of MAG, and in particular, it is necessary to clarify the 
role and place of the canonical energy-momentum tensor as a source of the 
gravitational field. Since one does not have the coframe (tetrad) among the 
fundamental variables, the corresponding field equation is absent. Here we 
demonstrate that one can always rearrange the field equations of MAG in 
such a way that the canonical energy-momentum tensor is recovered as one 
of the sources of the gravitational field. 

Assuming standard minimal coupling, the total Lagrangian density of 
interacting gravitational and matter fields reads 

-C = ^{9ij,Rijk\Nki^) + (51) 

In general, the gravitational Lagrangian density is constructed as a dif- 
feomorphism invariant function of the curvature, torsion, and nonmetricity. 
However, in view of the relations m and (EH), we can limit ourselves to 
Lagrangian functions that depend arbitrarily on the curvature and the dis¬ 
tortion tensors. The matter Lagrangian depends on the matter fields 
and their covariant derivatives (I34p . 

The field equations of metric-affine gravity can be written in several 
equivalent ways. The standard form is the set of the so-called “first” and 
“second” field equations. Using the modified covariant derivative defined by 
(EH), the field equations are given by 

^- eifc' = -V, (52) 

= 6'/. (53) 
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Here the generalized gravitational field momenta densities are introduced by 


:= 


2 


d‘^ 

dRkii^ ’ 






:= 


and the gravitational hypermomentum density 



= 


dNki^ ■ 


dQkij 


(54) 


(55) 


Furthermore, the generalized energy-momentum density of the gravitational 
field is 

Gifc' = 6i^ + ]^Qkin^'" + TkrSj'L + Rkln^^'^m. (56) 

The sources of the gravitational field are the canonical energy-momentum 
tensor and the canonical hypermomentum densities of matter, respectively: 


6 */ 


d2^raa.t ^ , A H O 


•^mat 


dT.J 






.B 


(57) 

(58) 


It is straightforward to verify that instead of the first field equation (I52p , 
one can use the so-called zeroth field equation which reads 


2 


^9ij 


t*i_ 


(59) 


On the right-hand side, the matter source is now represented by the metrical 
energy-momentum density which is defined by 


iij 2 


dg'-^ 


(60) 


The system ([52]) and ([53]l is completely equivalent to the system (l53]l and 
(|59p . and it is a matter of convenience which one to solve. 


2.3 Prom densities to tensors 

In actual metric-affine gravity models the Lagrangian densities are con¬ 
structed in terms of the Lagrangian functions: T = y/—gL, QJ = y/—gV, 
= 'J—gLma.t- Accordingly, we find for the gravitational field momenta 

(61) 

and j = ^^—gE^j . The gravitational sources (I57p , ([551) and ([Blip are rewrit¬ 
ten in terms of the canonical energy-momentum tensor, the canonical hyper¬ 
momentum tensor, and the metrical energy-momentum tensor, respectively: 

V = = (62) 
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The usual spin arises as the antisymmetric part of the hypermomentum, 

Ti/ := (63) 

whereas the trace = A*,^ is the dilation current. The symmetric traceless 
part describes the proper hypermomentum [25]. 

As a result, the metric-affine field equations (l5^ and are recast into 

+ = -T.k\ (64) 

+ = A*/. (65) 

2.3.1 Example: Matter model with microstructure 

In order to give an explicit example of physical matter with microstructure, 
we recall the hyperfluid model m- This is a direct generalization of the 
general relativistic ideal fluid variational theory |42l H3] and of the spinning 
fluid model of Weyssenhoff and Raabe |44| I45j . Using the variational princi¬ 
ple for the hyperfiuid [JT], one derives the canonical energy-momentum and 
hypermomentum tensors: 

Sfc' = EPk-p{6i-Evk), (66) 

= EJm^, (67) 

where u* is the 4-velocity of the fluid and p is the pressure. Fluid elements are 
characterized by their microstructural properties: the momentum density Pk 
and the intrinsic hypermomentum density Jm'^. 

3 General Lagrange-Noether analysis 

As a first step, we notice that the gravitational (geometrical) and material 
variables can be described together by means of a multiplet, which we denote 
by = {gij,Thi^ We do not specify the range of the multi-index J at 
this stage. The matter fields may include, besides the true material variables, 
also auxiliary fields such as Lagrange multipliers. With the help of the 
latter we can impose various constraints on the geometry of the spacetime. 
Furthermore, we can use the Lagrange multipliers to describe models in 
which the Lagrangian depends on arbitrary-order covariant derivatives of 
the curvature, torsion, and nonmetricity. Then the general action reads 

1 = j (68) 

where the Lagrangian density T = £(<h‘^, depends arbitrarily on the 

set of fields and their first derivatives. 

Our aim is to derive Noether identities that correspond to general coordi¬ 
nate transformations. However, it is more convenient to start with arbitrary 
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infinitesimal transformations of the spacetime coordinates and the matter 
fields. They are given as follows: 


x'\x) = X® + fe®, 




(69) 

(70) 


Within the present context it is not important whether this is a symmetry 
transformation under the action of any specific group. The total variation 
dZOD is a result of the change of the form of the functions and of the change 
induced by the transformation of the spacetime coordinates ([6^ . To distin¬ 
guish the two pieces in the field transformation, it is convenient to introduce 
the substantial variation: 




(71) 


By definition, the substantial variation commutes with the partial derivative, 
We need the total variation of the action: 


61 = j [d^x 62 + 6{d*x) £] . 


(72) 


A standard derivation shows that under the action of the transformation 
(f69]) - (f70]l the total variation reads 


61 = / d^x 


62 


82 






Here the variational derivative is defined, as usual, by 




82 


6^j ■ 


- 8i 


82 


8{8i^-^) 


(73) 


(74) 


3.1 General coordinate invariance 


Now we specialize to general coordinate transformations. For infinitesimal 
changes of the spacetime coordinates and (matter and gravity) fields (l6^ 
and dZOD we have x® —x® -|- (5x®, gij —> gij + 6gij, T/jjl —> T^jl -|- (iT^jl, 
and =>■ where variations are given by ([2]), ([3]), ([4]), and ([32]). 

Substituting these variations into (I73]l . and making use of the substantial 
derivative definition (|71l) . we find 


61 



e + {8ie) + {dfje) 


(75) 
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where explicitly 
nk = 


^9ij 
+ di 


dk9ij “1“ 


(5£ ^ 82, ^ 

+ TT;:;— —dk^ln + 


82 


Spj 


8Tin- 


88iTin^ 


8k8iTir, 


‘‘k — 


^ 88igrnn 

^ 52 52 

- 9k j + 


82 82 j , 

8k9mn + T^HTTr ^k'P “ 4-^ > 


88i'p^ 


{(y^B)k'P^ + 


( 76 ) 


^9ij 
82 


5il;^ 

8k9mn “1“ 8j 


82 

88i'ip^ 
9nk + 


dki^^ - Si2 


^82 92 i Q 

2— - Qnk + ■7T7r—ir{cr^B)k W 


88igmn"88j gin'’' 5 ^ 4 ^ 

- Tik^+^ Vkp - T;/+ 


82 


8Tu^ 

82 


Mi- il 


8Ti8 

■ 8n2kl^ + 


8T 


82 


= 


82 


88^i'P^ 

82 

+ 2 




88nTu'^ 

82 


ij 

■ - 


In 
82 

88nTlm^ 


■ 8k2lr, 


8nTlm\ (77) 


+ 


82 


kl 




98(^igj'jn 

82 


9%r,.)r 

82 „ ^ 82 

Ik c\c\ T-i k ^ 


9kn + r,,'" - 


A8 \l\j) 


dd(^ir\in\' 


d8^n^ij)’^' 


(78) 

(79) 


If the action is invariant under general coordinate transformations, 51 = 0, 
in view of the arbitrariness of the function and its derivatives, we find the 
set of four Noether identities: 


= 0, fifc* = 0, nk^ = 0, Qk^n ^ 0. (80) 


General coordinate invariance is a natural consequence of the fact that 
the action (|68p and the Lagrangian 2 are constructed only from covariant 
objects. Namely, £ = 2{'ip^,V, gij, Rkii^, Nk/) is a function of the met¬ 
ric, the curvature ([5]), the torsion (l6|), the matter fields, and their covariant 
derivatives (fM)l . Denoting 


:= 


82 


8R 


ijk 


I ’ 


:= 


82 

8Nij^ ’ 


we find for the derivatives of the Lagrangian 


82 


82 

sajy 

82 

88kgij 


+2p^^^krnl^ +2p^9^Ynk\ 

2p9\, 


(81) 


(82) 

(83) 

(84) 


As a result, we straightforwardly verify that = 0 and = 0 are 

indeed satisfied identically. 
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Using (|82l) and (|83D . we then recast the two remaining Noether identities 
([76]) and dZZ]) into 


+ V, (cT^bU^ Tkn^ + P^^^mdkRiln^ 

+ Rlkn^ + 

+ dkOmn, (85) 

I o D m I ^Ini D m Anm d i 

+ m-ti^kln “T p m^lnk P k-^lnm 

+ Vn 5ifc = 0. (86) 

Notice that the variational derivative (I74p with respect to the matter 
fields can be identically rewritten as 


8^) • 

and turns out to be a covariant tensor density. It is also worthwhile to 
note, that the variational derivative w.r.t. the metric is explicitly a covariant 
density. This follows from the fact that the Lagrangian depends on gij not 
only directly, but also through the objects Qkij and Taking this into 

account, we find 


(IT _ ^ / 

6ip^ ^ \ 


(5£ 

^9ij 


j On 

dgij 

dS: 1 - 

-Vr 

dg^, 2 


d£ 

ddngij 

g{ni)j 


+ 


( 88 ) 


The Noether identity (|86p is a covariant relation. In contrast, (|85p 
is apparently non-covariant. However, this can be easily repaired by re¬ 
placing Hfc = 0 with an equivalent covariant Noether identity: Q'k ■= 
— ^kn^^rri^ = 0. Explicitly, we find 




(5£ 


Vktp^ + Vi 


dH 


dVi^jj^ 


Vk^^ - 6l£ 


g ST ^B 


+ 


dViipA 

62. 1 


dViil^A 


_ iv _|_ ,Xnj)i _ ..{ij)n\ 

6g,j ^ ) 


Qk 




+ P^^mVkNin^ + P^,r^VkRiln^ 


(89) 
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On-shell, i.e. assuming that the matter fields satisfy the field equations 
6£,/6'ip^ = 0, the Noether identities ([86]) and (1891) reduce to the conservation 
laws for the energy-momentum and hypermomentum, respectively. 


Equation (1861) contains a relation between the canonical and the metrical 


energy-momentum tensor, and the conservation law of the hypermomentum. 
In the next section we turn to the discussion of models with general non- 
minimal coupling. 

4 Conservation laws in models with nonminimal 
coupling 

The results obtained in the previous section are applicable to any theory in 
which the Lagrangian depends arbitrarily on the matter fields and the grav¬ 
itational field strengths. Now we specialize to the class of models described 
by an interaction Lagrangian of the form 


L = F{gij,Rki^,Nkn V.V’^). 


(90) 


Here is the ordinary matter Lagrangian. We call F = 

F{gij,Rkii^,Nki') the coupling function and assume that it can depend ar¬ 
bitrarily on its arguments, i.e., on all covariant gravitational field variables 
of MAG. When E = 1 we recover the minimal coupling case. 

4.1 Identities for the nonminimal coupling function 

As a preliminary step, let us derive identities which are satisfied for the 
nonminimal coupling function F = F{gij, Ri^n^, For this, we apply 

the above Lagrange-Noether machinery to the auxiliary Lagrangian density 
-^0 = V~9 This quantity does not depend on the matter fields, and both 
(|86|) and (I89|) are considerably simplified. In particular, we have 



(91) 


Then we immediately see that ([5U|) and (15^ reduce to 



(92) 





(93) 


Here we denoted 



( 94 ) 
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Notice that for any tensor density we have relations (I30|) and (1291) . 

The identity (j92p is naturally interpreted as a generally covariant gen¬ 
eralization of the chain rule for the total derivative of a function of several 
variables. This becomes obvious when we notice that (1881) implies 


_ _v„ 

2 




dF 


Qkij — o k 9 ij ■ ( 95 ) 


It should be stressed that (1^21) and ([USP are true identities, they are satisfied 
for any function F{gij, Rkii , N'fcz*) irrespectively of the field equations that 
can be derived from the corresponding action. 


4.2 Conservation laws 


Now we are in a position to derive the conservation laws for the general 
nonminimal coupling model (1901) . Recall the definitions of the canonical 
energy-momentum tensor (j57p . the canonical hypermomentum tensor (|58l) 
and the metrical energy-momentum tensor (I60p . 

In view of the product structure of the Lagrangian (fTOp . the derivatives 
are easily evaluated, and the conservation laws (fHHp and (fMP reduce to 


- F4* - Vn (TA'fc-) + FSfc* + 


2Fk 


I n^iln 75 m I ^Ini 75 m ^Inm 75 i 

T m-^kln T P m^lnk P k-^lnm 




^li 


FF^kNin^ -P^\Nik 

- Vn ( I Qjk 


Tmat — 0, 


1 


Vi (FSfc') + F + -f^Qk 




+ {p^^^m^kRiln^ + P^^m^kNln'^ + 

Q,. .| 


— pi] 


Tmat — 0. 


(96) 


(97) 


After we take into account the identities ((921) and (f93p . the conservation 
laws ((Mp and (1^ are brought to the final form: 


FSfc* = W + V^FAV), (98) 

Vi (FSfc') = F - A^jRkim" - - TmatV^T. (99) 

Lowering the index in (1981) and antisymmetrizing, we derive the conservation 
law for the spin 

-|- Vn {FTij^) + Qnl[iA^= 0. (100) 

This is a generalization of the usual conservation law of the total angular 
momentum for the case of nonminimal coupling. 
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4.3 Rewriting the conservation laws 

Using the definition (I3ip and decomposing the connection into the Rieman- 
nian and non-Riemannian parts (jlOp . we can recast the conservation law 
(j98p into an equivalent form 

- 4 ' + ( 101 ) 

In a similar way we can rewrite the conservation law (j99p . At first, with the 
help of fjl2p and m we notice that 

F = F{ti^ - J:nNk/ + FJ^i^Nik'. (102) 

Then substituting here (llOip and making use of (fTop , (fMIl , and the curvature 
decomposition (fTTP . after some algebra we recast ([99P into 

V,{F(Sfc^+AViVfcm”)} = -FA^Jih^m"-^kNimn-Lm.tVkF (103) 

For the minimal coupling case, such a conservation law was derived in |46l 
ST]. The importance of this form of the energy-momentum conservation 
law lies in the clear separation of the Riemannian and non-Riemannian 
geometrical variables. As we see, the post-Riemannian geometry enters (jl03l) 
only in the form of the distortion tensor Afcj* which is coupled only to the 
hypermomentum current This means that, in the minimal coupling 

case, ordinary matter - i.e. without microstructure, A”^„® = 0 - does not 
couple to the non-Riemannian geometry. In contrast, in the nonminimal 
coupling case, the derivative of the coupling coupling function F on the 
right-hand side of (jl03p may lead to a coupling between non-Riemannian 
structures and ordinary matter. 


4.4 Conserved current induced by a spacetime symmetry 

Every generalized Killing vector field generates a conserved current. This 
can be demonstrated from the analysis of the system (llOip and (jl03p as 
follows. Let us contract equation (IIOII) with V and equation (11031) 
with , then subtract the resulting expressions. Note that the contrac¬ 
tion = 0 vanishes because the first factor is a symmetric tensor and 

the second one is skew-symmetric. Then after some algebra we derive 


V, = FA^JC^N^m^ - 


Here we associate a current with a Killing vector field via 


/* := F 


^ m I 
n 


(104) 


(105) 


Note that the transposed connection appears here. The right-hand side 
of (jl04p depends linearly on the Lie derivatives along the Killing vector: 
C^F = C^VkF and 

CcNim" = + {ViC^)Nkm" + {VmC’^)Nik^ - {VkC)Nim^. (106) 
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When is a generalized Killing vector, we have = 0 in view of 

(I46p . Furthermore, recalling that F = F{gij, ^ Nj^j^)^ we find 

by making use of (I46|) . (HTll . and (jiQl) . 

As a result, the right-hand side of (I1U4D vanishes for the generalized 
Killing vector field, and we conclude that the induced current (|105p is con¬ 
served 

vj* = 0. (108) 

This generalizes the earlier results reported in [3011481H9] . In section [6.1.2l we 
will show that there is a conserved quantity constructed from the multipole 
moments which is a direct counterpart of the induced current (|105p . It is 
worthwhile to give an equivalent form of the latter: 


/* = F 


- {VmCW 


(109) 


4.5 Riemannian limit 

Our results contain the Riemannian theory as a special case. Suppose that 
the torsion and the nonmetricity are absent Tij^ = 0, Qkij = 0, hence Nij^ = 
0. Then for usual matter without microstructure (i.e. matter with = 0) 

the canonical and the metrical energy-momentum tensors coincide, = 
tk- As a result, the conservation law (IMp reduces to 

VRfc® = i (-Lmat^i - 4') ViF. (110) 

This conservation law for the general nonminimal coupling model was de¬ 
rived earlier in [26] without using the Noether theorem, directly from the 
field equation^. The old result established the conservation law for the case 
in which F = F{Rijk ) depends arbitrarily on the components of the curva¬ 
ture tensor, correcting some erroneous derivations in the literature, see [26] 
for details. 

Quite remarkably, (IllOl) generalizes the earlier result to the case in which 
the nonminimal coupling function F is a general scalar function of the cur¬ 
vature tensor. 


5 General multipolar framework 

In this section we derive “master equations of motion” for a general extended 
test body, which is characterized by a set of currents 




( 111 ) 


®Notice a different conventional sign, as compared to our previous work | 26 | . 
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Normally, these are so-called Noether currents that correspond to an invari¬ 
ance of the action under a certain symmetry group. However, this is not 
necessary, and any set of currents is formally allowed. We call dynam¬ 
ical currents. The generalized index (capital Latin letters A,B,...) labels 
different components of the currents. 

As the starting point for the derivation of the equations of motion for 
generalized multipole moments, we consider the following conservation law: 

( 112 ) 

On the right-hand side, we introduce objects that can be called material 
currents 

(113) 

to distinguish them from the dynamical currents The number of com¬ 
ponents of the dynamical and material currents is different, hence we use a 
different index with dot the range of which does not coincide with 

that oi A, B,. At this stage we do not specify the ranges of both types 
of indices, this will be done for the particular examples which we analyze 
later. As usual, Einstein’s summation rule over repeated indices is assumed 
for the generalized indices as well as for coordinate indices. 

Both sets of currents and are constructed from the variables 
that describe the structure and the properties of matter inside the body. In 
contrast, the objects 

(114) 

do not depend on the matter, but they are functions of the external classical 
fields which act on the body and thereby determine its motion. The list of 
such external fields include the electromagnetic, gravitational and scalar 
fields. 

We will now derive the equations of motion of a test body by utilizing 
the covariant expansion method of Synge [7j . For this we need the following 
auxiliary formula for the absolute derivative of the integral of an arbitrary 
bitensor density (x, y) (the latter is a tensorial function of 

two spacetime points): 

£ j = j j {lib) 

S(s) T,{s) E(s) 


Here := dxy^/ds, s is the proper time, ^ and the integral is 

performed over a spatial hypersurface. Note that in our notation the point 
to which the index of a bitensor belongs can be directly read from the index 
itself; e.g., yn denotes indices at the point y. Furthermore, we will now 
associate the point y with the world-line of the test body under consider¬ 
ation. Here a denotes Synge’s [7] world-function and its first covariant 
derivative; gyx is the parallel propagator for vectors. For objects with more 
complicated tensorial properties the parallel propagator is straightforwardly 
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generalized to x and G^ We will need these generalized propagators 
to deal with the dynamical and material currents and K^. More details 
are collected in the appendix on our conventions. 

After these preliminaries, we introduce integrated moments for the two 
types of currents via (for n = 0,1,...) 

jj/i-ynib = (116) 

E(r) 

^yi-yr^Yoy' ^ {-lYjay^---ay-G^°Xo9^'x'2^°'"'w^''(n:^», (117) 

E(r) 

E(r) 


Here the densities are constructed from the currents: 

and = y/—gK^°. Integrating (jll2p and making use of (I115j) . we find 
the following “master equation of motion” for the generalized multipole 
moments: 


D 

ds 


jyi-y-nYo _ _^y(yijy2-yn)Yo ^j(jji...yn-i\Yo\yu) 

-7\'y"yr.+, {iy^-y-^'y" +jy^-y-^'vy' 
-A,,>Y"^oiy^-y-^''y' - K>Y"^°.„„^Ay^-y^+^^''y‘ 


y'Y" \yn+Y 

^ lyri+i 


Vj/'y" 

V 


+En 


k=2 


k\ 


Vn + l • • -yn+k ^ 


iy2—yn)yn+i---yn+kYoy' 




_A ,,Yo jyi-yn+kY''y' 

l\.ytY" ■,yn+l...yn+k'' 


-n^o > 


''iyn + l'"yn+k 


m 


yi'"yn+k^ 


(119) 


5.1 Electrodynamics in Minkowski spacetime 


To see how the general formalism works, let us consider the motion of elec¬ 
trically charged extended test bodies under the influence of the electromag¬ 
netic field in flat Minkowski spacetime. This problem was earlier analyzed 
by means of a different approach in |50j . 

In this case, it is convenient to recast the set of dynamical currents into 
the form of a column 


= 


j^kj 


( 120 ) 
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where is the electric current and is the energy-momentum tensor. 
Physically, the structure of the dynamical current is crystal clear: the mat¬ 
ter elements of an extended body are characterized by the two types of 
“charges”, the electrical charge (the upper component) and the mass (the 
lower component). 

The generalized conservation law comprises two components of different 
tensor dimension: 

^3 ( yI^3 ) ^ ( 

where the lower component of the right-hand side describes the usual Lorentz 
force. 

Accordingly, we indeed recover for the dynamical current (jl20h the con¬ 
servation law in the form ()112l) . where = 0 and 



( 122 ) 


The generalized moments (I116l) - (|118l) have the same column structure, re¬ 
flecting the two physical charges of matter: 


, / jyi---yny' \ 

^ j , ,124) 


whereas = 0. 

As a result, the master equation (jllQIl reduces to the coupled system of 
the two sets of equations for the moments: 


I^jVvyn _ _^y{yijy2--yn) ^j^{yi---yn)^ 

ds"' ’ 


D 

ds 


pyvynyo _ _j^y(yipy2---yn)yo y^j^(yi-yn-i\yo\yn) 


_F, 


,yo 


k=l 


y yn+l—yn+k 


(125) 


^yi-yn+ky ( 126 ) 


These equations should be compared to those of [50] 


6 Equations of motion in metric-affine gravity 

We are now in a position to derive the equations of motion for extended 
test bodies in metric-affine gravity. As a preliminary step, we rewrite the 
conservation laws ([Mj) and ([Mil as [22] : 



TT m Am J . 1 . i 

^jm n ~r ^k ^k ? 

(127) 


'{/ri^ j R n Am j ^ D 3 AT 

k‘^n ^kjm ^ n ^^^kj -^/c-^mat* 

(128) 
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Here we denoted := log-F, and 


Vjn^ := H,,5^ + 7VV- 

Introducing the dynamical current 



and the material current 



(129) 

(130) 


(131) 


(132) 


we then recast the system (I127p and (jl28|) into the generic conservation law 
(11121) . where we now have 

AjB 

H"^ ■ 

Like in the previous example of an electrically charged body, the matter 
elements in metric-affine gravity are also characterized by two “charges”: the 
canonical hypermomentum (upper component) and the canonical energy- 
momentum (lower component). This is reflected in the column structure of 
the dynamical current (I13ip . The material current (I132p takes into account 
the metrical energy-momentum and the matter Lagrangian related to the 
nonminimal coupling. The multi-index A = {ik,k}, whereas A = {ik,l}- 
Accordingly, the generalized propagator reads 



-5)5% 

ji'k' 





^4^ ) ■ 


(133) 

(134) 




nVi ny2 

y xiU X 2 

0 

0 

9^^ XI 


(135) 


and we easily construct the expansion coefficients of its derivatives from the 
corresponding expansions of the derivatives of vector propagator 


vT), 


where we denoted 


- 


0 

yk+2 

V 0 

'y^°y'y2---yk+2 


^{yoy} . , — ^yo , xy ^ ^y ,, xyo 

^ {y y }y2---yk+2 ~ ' y'y2-yk+2^y" ^ / y”y2-yk+2^y' ■ 

In particular, for the hrst expansion coefficient (k = 1), we hnd 


'y^^°^\y'y"}y2y3 


y'y2y3^y" + ^^y''y2y3^y° 


^yo , — -pyo , 

/ yy2y3 ~ 2 yy2y3- 


(136) 


(137) 


(138) 

(139) 
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For completeness, let us write down also another generalized propagator 


■ - 
X - 


y XiU X2 

0 

0 

1 


(140) 


The last step is to write the generalized moments (|116p - (lll8l) in terms 
of their components: 


jVi-ynY 

jyi-'-VnYyo 

jj^yi-'-ynY 


hyvyny'y” 


P' 


yi---yny 


qy 

j^yi-yny'yo 




^yi--- 


Vn 


(141) 

(142) 

(143) 


For the two most important moments, “h” stands for the hypermomentum, 
whereas “p” stands for the momentumJl Finally, substituting all the above 
into the “master equation” (11191) . we obtain the system of multipolar equa¬ 
tions of motion for extended test bodies in metric-affine gravity: 


_ uyi-ynyayb 

ds 


—n y(y^h^'2-"y^')y<^yb _(- ^ q(yi---yn-i\yayb\yn) 
_i^j^yi---ynybya _ ^yi—ynyayb 

f^qyi-yn+iy’yby" ^ yy\yi-yn+iy'yb 

1 


—TT , „yaybnyi---yny y yo 

0 ir\'7r t r ' 


-_Ryby,y„^^^^f^qyi-yn+iyayy ^ ^y j^yi-yn+iyay 

y yu_ TT , ..yaVb ^yi-'-yn+iy y yo 

y ^yoy'y" lyn+iH 


'yoy y 

oo 




k=2 


k\ 


fl(yi , hy2-yn)yn+l-yn+kyayb 

V r y y-n+l-'-ynA-k' 


-t-( 1) l^°'y'y"y„+i...y„+k 
+ (“1) l^’’y'y"yn+i---yn+k 


qyi-'-yn+ky'yby" _|_ yy”Ryi---yn+ky'yb'^ 
qyi-yn+kyay'y” _|_ yy”Ryi---yn+kyay''^ 


— 1—/ gy2-yu)yn+i.-y„+kyayby' 

V y ^71 + 1'"^/n+fc 


—TT , ,,y°-yb „yi---yn+ky'y''yo 

'^yoyy iyu+i-.-yn+k^T 


(144) 


^Note that in order to facilitate the comparison with our previous work [sa, we provide 
in appendix the explicit form of integrated conservation laws (fml) and dmj, as well as 
the generalized integrated moments HU} - diia} in the notation used in m- 
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^ „yi-ynya 

ds^ 


_^^(yipj/2---j/n)ya _|_ j^j^{yi-yn-i\ya\yn) 

_ j\ya^yi---yn _ A^°--y 

— V„ ,y<^ky^---yny'y'' —v„ ,y<^ uyi--yn+iy'y" 

Vyllyl K Vyllyl 

--KVa , ,, (uyi-'-yn+iy'y" \ „,y''yi-yn+iy' 

2^ yy yn+i\^ ^ 

— 7??^“ , „ay^---yny'y”yo _ 

yoy y 'i 

— \nyo. „ , ..yi-yny'y" _ 1/0^“ „ , ,,yi'" 2 /n+i?;'?;" 

U/ y'/y'fj^ ^ y"y>-,yn+lh‘‘ 


, ,, riyi-"yn+iy'y"yo 

yoy y ,yn+i'i 


+e4 


k=2 


k\ 




-\_(^ A^ycL , „ (uyi---yn+ky'y'' \ .,y''r,y^---y^+i^y'\ 

-*-7 7 y'y"yn+i---yn+k -r u p j 


y y yn+i---yn+k 
k - {yi 


-{-Ifna^y^y: 


-Ry^ 


y'yn+i.-yn+k 


yoy'y”-,y 71+1 ■■■yn+k^ 


j^y2---yn)yn+i---yn+kyay' 


_T/,, ,yo. uyi---yn+ky y 

^y y Wn+i-.-yn+k'^ 


_ ny<^ „ , ..yi-yn+ky y 

2 ^ y y ,yn+i---yn+kt^ 


_ AVa cyi-'-Vn+k 

?2/n + l •••^n+fc 


(145) 


6.1 Special cases 

The general equations of motion (11441) and (|145p are valid to any multipolar 
order. In the following sections we focus on some special cases, in particular 
we work out the two lowest multipolar orders of approximation, and consider 
the explicit form of the equations of motion in special geometries. 

6.1.1 General pole-dipole equations of motion 

From (I144p and (I145p . we can derive the general pole-dipole equations of 
motion. The relevant moments to be kept at this order of approximation 
are: A:“^, and Since all objects are now 

evaluated on the world-line, we switch back to the usual tensor notation. 
For n = 1 and n = 0, eq. (I144p yields 

0 = _ ^abc ^^bca _^a^bc^ 

Dj^ab ^ J^ba _ ^ab _jj^^^ab^dec_ 
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Furthermore for n = 2,1,0 equation (|145p yields 


0 

_ j^ia\c\b) _ y{a.pb)c^ 

(148) 

as 

= A:'“ - uV - 

(149) 


= - R^dbcq^^'^ - \Q\b^^^'^ 



-A% - ]^R\db 

T/ a ubcd ^/na ,.hcd a a 

^dc 2^^ dc\bl^ ;&S • 

(150) 


Rewriting the equations of motion Let us decompose (11461) and (jl47p 
into symmetric and skew-symmetric parts: 


pabc ^ 

J^a{bc) ^ ^{bc)a _ ^aj^(bc) ^ 

(151) 

0 = 


(152) 

= 

ds 

(153) 

P_h['^b\ ^ 

ds 


(154) 


As a result, we can express the moments symmetric in the two last indices 
^ab _ ^{ab) ^cab _ ^c{ab) general, this is possible also for an arbitrary 
order ^^<=1 = ^ci...cn(ab}^ gf i^j^g other moments. 

Let us denote the skew-symmetric part := as this simplifies 

greatly the subsequent manipulations and the comparison with m- 

The system of the two equations (|148p and ()152l) can be resolved in terms 
of the 3rd rank /c-moment. The result reads explicitly 


j^abc 


_j_^[afe]c ^[ac]b q[bc]a^ 



This yields some useful relations: 

J^albc] ^ _ya^bc^^[bc]a^ 

J.lab]c ^ ^ f^plab] _ ^abj ^ ^[a6]c_ 


The next step is to use the equations (I15ip . (jl53h together with (jl55h 
and to substitute the /r-moments and A;-moments into (11471) and (|149l) - (jl50p . 
This yields a system that depends only on the p, h, q, and ^ moments. 

Let us start with the analysis of (|150p . The latter contains the combina¬ 
tion + y[<ipb]c -^iigpg the skew symmetry is imposed by the contraction 

with the Riemann curvature tensor, which is antisymmetric in the two last 
indices. Making use of (|155p . we derive 


+ y[bpa]c ^ ^abc ^ ^acb _ ^bca^ 


(158) 
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where we introduced the abbreviation 


:= v'^ ( 


(159) 


Note that by construction 

Then by making use of the Ricci identity we find 


- \R\db ( 


= 


bed 


^[cd]b ^ ^b / [cd] _ ^cd 


(160) 


Substituting from (I149p and from (|153p . we find after some algebra 

- - (p“ + v^A, 

-yl“ Ab - + {N\bNdc" - N\nNd^b) A- ( 161 ) 


Further simplification is achieved by noticing that 

DA DA 

v^’Ab = —, k’^^‘^A,,A, = (162) 

as as 

where we used (|148p and recalled that A^ = A-ij. 

Analogously, taking k^^^^ from (I156p and from (jl51|) . we derive 

- = - A,,X^^ + A - (163) 

We can again use A), = A-b and (jl48p to simplify 

- A,.,A = (164) 

as 

After these preliminary calculations, we substitute (I160I) - (I164I1 into (I150p 
to recast the latter into 


- (Fp“ + FN\dh^+p’^^VbF^ = FR\,dV^ - 5^“^) 


—Fq' 


-FA 


cbd 


R^dcb - R^dcb - N\b-,d - N\bNd" + N^cnN/^b 




(165) 


Finally, combining (jl47p and (I149h to eliminate k^°' we derive the equa¬ 
tion 

£('^afe_^ab^ = A-^“(/ + AA) 

+ A Ad - AlVdc“ + A Ac 
-r7lA(g“^^-0^c. (166) 


Following m, we introduce the total orbital and the total spin angular 
moments 

L“^=2p[“''l, A :=-2/i[“''l, (167) 
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and define the generalized total energy-momentum 4-vector and the gener¬ 
alized total angular momentum by 



:= F(p“ + A%rfh^'^) + 

(168) 


77“^ := F(L“^-F5“^). 

(169) 

Then, taking into account the identity (I14p. which with the help 
raising and lowering of indices can be recast into 

of the 


= -R^dcb + R'^dcb + N\b-,d + N\bNdc^ - Nd^bN\n, 

(170) 

we rewrite the pole-dipole equations of motion (|165p and (|166p in the final 

form 



DV^ 

ds 

= \R\cdv'^J^ + Fq^^^V^^Ndcb 


ds 


(171) 



_2^[ay6]^_ 

(172) 


The last equation arises as the skew-symmetric part of (|166ll . whereas the 
symmetric part of the latter is a non-dynamical relation that determines the 
^ab j^ioment 

1 / N 

^ __ ^ _^(a j-pfe) ^ ^ ^{a^b) 

+ (gHfe + qMa _ q{^b)c^A^. (173) 

Here the symmetric moment of the total hypermomentum is introduced via 

^ab ^{ab) _ j^iab) _ 


6.1.2 Conserved quantity 

The equations of motion for the multipole moments are derived from the 
conservation laws of the energy-momentum and the hypermomentum cur¬ 
rents Sfc* and A”*„*. In section 14.41 we demonstrated that every generalized 
Killing vector induces a conserved current constructed from Sfc® and 
Quite remarkably, there is a direct counterpart of such an induced current 
built from the multipole moments. 

Let C,^ be a generalized Killing vector, and let us contract equation (11711) 
with Qa and equation (|172p with aCb^ and then take the sum. This yields 

£ (^“Ca + = Fq^^^C^Ndcb e^aFcF (175) 
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On the right-hand side, the Lie derivatives of the distortion tensor and of 
the coupling function both vanish in view of (|16]) and (I107p . 

Consequently, we conclude that for every generalized Killing vector field 
the quantity 

V\a + = const (176) 

is conserved along the trajectory of an extended body. 

We thus observe a complete consistency between (I104p . (llOSp . (I109p and 

dlTel), (UZH]). 

6.1.3 Coupling to the post-Riemannian geometry: Fine structure 

Let us look more carefully at how the post-Riemannian pieces of the grav¬ 
itational field couple to extended test bodies. At first, we notice that the 
generalized energy-momentum vector (|168l] contains the term that 

describes the direct interaction of the distortion (torsion plus nonmetricity) 
with the intrinsic dipole moment of the hypermomentum. Decomposing 
the latter into the skew-symmetric (spin) part and the symmetric (proper 
hypermomentum -|- dilation) part, we find 

(177) 

Here we made use of (jl3p . This is well consistent with the gauge-theoretic 
structure of metric-affine gravity. The second term shows that the intrinsic 
proper hypermomentum and the dilation moment couple to the nonmetric¬ 
ity, whereas the first term displays the typical spin-torsion coupling. 


Similar observations can be made for the coupling of higher moments 
which aouear on the rieht-hand sides of 1I171D and (I172D - and thus determine 
the force and torque acting on an extended body due to the post-Riemannian 
gravitational field. In order to see this, let us introduce the decomposition 

^ _|_ qcab 

(178) 

into the two pieces 


qO-bc 1 ^glac]b ^[bc]a _ ^{ab)c^ ^ 

(179) 

qabc 1 ^qlab]c ^[ac]b _ ^[fec]a^ _ 

(180) 


The overscript “d” and “s” notation shows the relevance of these objects to 
the dilation plus proper hypermomentum and to the spin, respectively. By 
construction, we have the following algebraic properties 


q[o-b\c = 0, = 0. (181) 
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Making use of the decomposition (11781) and of the explicit structure of 
the distortion m, we then recast the equations of motion (I17ip and (jl72l) 
into 

- eV“F - ^ (182) 

-2v 

_2^[“V^1f. (183) 

Now we clearly see the fine structure of the coupling of extended bodies to 
the post-Riemannian geometry. The first lines in the equations of motion 
describe the usual Mathisson-Papapetrou force and torque. They depend on 
the Riemannian geometry only. A body with the nontrivial moment (jlSOj) 
is affected by the torsion field, whereas the nontrivial moment (|179l) feels 
the nonmetricity. This explains the different physical meaning of the higher 
moments (jl79p and (jl80l) . In addition, the last lines in (I182p and (I183p 
describe contributions due to the nonminimal coupling. 



ds 


ds 


6.1.4 General monopolar equations of motion 

At the monopolar order we have nontrivial moments p“, and The 

nontrivial equations of motion then arise from (11440 for n = 0 and from 
(jl450 for n = 1, n = 0: 


0 


(184) 

0 


(185) 

ds 


(186) 

The first two equations 

(I184P and (I185p vield 



= 0, uiy] = 0, 

(187) 

and substituting (jl84P. 

(11851) and (|187p into (11861) we find 



D(Fp°‘) 

ds 

(188) 


From (I187p we have with the mass M := v°‘Pa, and this allows to 

recast (jl88p into the final form 


M 




ds 


= -e(r-^V)- 


F 


(189) 
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Hence, in general the motion of nonminimally coupled monopole test bod¬ 
ies is nongeodetic. Furthermore, the general monopole equation of mo¬ 
tion (jl89l) reveals an interesting feature of theories with nonminimal cou¬ 
pling. There is an “indirect” coupling, i.e. through the coupling function 
F{gij, Rijk\Tij^, Qij^), of post-Riemannian spacetime features to structure¬ 
less test bodies. 

6.1.5 Weyl-Cartan spacetime 

In Weyl-Cartan spacetime the nonmetricity reads Qkij = QkSij^ where Qk 
is the Weyl covector. Hence the distortion is given by 

Nki = Kki + ^ (Q%J - Qk5) - QjSl) . (190) 

The contortion tensor is constructed from the torsion, 

Kk/ = - ^ {Tki + T\j + Tjk)- (191) 

As a result, the generalized momentum (I168p in Weyl-Cartan spacetime 
takes the form: 

= Fp^-^ (K^cdS'^^ - QbS^^ + (192) 

Here we introduced the intrinsic dilation moment D := gabh°'^- 

Substituting the distortion ()190p into ()17ip and ()172p . we find the pole- 


dipole equations of motion in Weyl-Cartan spacetime: 


DV^ 

ds 

= + Fq^^^V^T,bd 



+ Z^V^Qb - - C^VbV^F, 

(193) 

Djab 

ds 



+2FZ^^Q^^ - 

(194) 

Here we introduced the trace of the modified moment (I179p 



_ „ fibca _ ^ „a,bc\ 

■— 9bc9 — ^Qbc (^9 - q - q J ■ 

(195) 


It is coupled to the Weyl nonmetricity. 


6.1.6 Weyl spacetime 

Weyl spacetime [52] is obtained as a special case of the results above for 
vanishing torsion. Hence the contortion is trivial 

Kabc = 0. (196) 

Taking this into account, the generalized momentum (I192p and the equations 
of motion (|193p and (jl94p are simplified even further. 

It is interesting to note that besides a direct coupling of the dilation 
moment to the Weyl nonmetricity on the right-hand sides of (|193h and (I194p , 
there is also a nontrivial coupling of the spin to the nonmetricity in (11921) . 
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6.1.7 Riemann-Cartan spacetime 

Another special case is obtained when the Weyl vector vanishes Qa = 0. 
Equations ()192p - (jl94l) then reproduce in a covariant way the findings of 
Yasskin and Stoeger |53] when the coupling is minimal {F = 1). For non- 
minimal coupling we recover our earlier results in [51j . 

7 Equations of motion in scalar-tensor theories 

The geometrical arena of scalar-tensor theories is the Riemannian space- 
time, hence = 0 which means that both the torsion Tij^ = 0 and the 
nonmetricity Qkij = 0 vanish. 

Scalar-tensor theories have a long history and they belong to the most 
straightforward generalizations of Einstein’s general relativity theory. In the 
so-called Brans-Dicke theory [SH [55l [56l EH [58] a scalar field is introduced 
as a variable “gravitational coupling constant” (which is thus more correctly 
called a “gravitational coupling function”). Similar formalisms were devel¬ 
oped earlier by Jordan [591 I60j . Thiry |61] and their collaborators using 
the 5-dimensional Kaluza-Klein approach. An overview of the history and 
developments of scalar-tensor theories can be found in (62] |63l [Ml ES] • 

Surprisingly little attention was paid to the equations of motion of ex¬ 
tended test bodies in scalar-tensor theories. Some early discussions can be 
found in [56l [Ml ET] , and in [68] the dynamics of compact bodies was thor¬ 
oughly studied in the framework of the post-Newtonian formalism. 

7.1 Scalar-tensor theory: Jordan frame 

We consider the class of scalar-tensor theories along the lines of |68j where 

the action I = f d'^x £ is constructed on the manifold with the spacetime 

j J J J 

metric dij. The Lagrangian density £ = y/—gL has the following form 

-2U^ + (197) 

This action is an extension of standard Brans-Dicke theory [M] to the case 
when we have a multiplet of scalar fields (capital Latin indices A,B,C = 
1,... ,N label the components of the multiplet). Here k = STrGjd^ denotes 
Einstein’s gravitational constant and in general we have several functions of 
scalar fields, 


J 1 / ^ J- J A 

L = ^(3) + g^^lAsdi^^dj 


F = F{^^), u = U{g^^), 1 ab = 1ab{^'^)- (198) 

J 

The Lagrangian dip, gij) depends on the matter helds xp and the grav¬ 

itational field. 
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The metric dij determines angles and intervals in the Jordan reference 

~ J 

frame. The Riemannian curvature scalar R{9) is constructed from the Jor¬ 
dan metric. With the help of the conformal transformation 


ij ^ 9ij — R 9ij 


(199) 


we obtain a different metric on the spacetime manifold. This is called an 
Einstein reference frame. 


7.2 Conservation laws: Einstein frame 

In the Einstein reference frame the Lagrangian density in the scalar-tensor 
theory reads £ = y/^L with 


1 


B 


L =— [-R +dj(p - 217) -t E gij). (200) 


1 


p-2. 




Here the scalar curvature R{g) is constructed from the Einstein metric gij, 
and 

7AB = -^{Iab + ^F,aF,b), U = -^u . (201) 

The metrical energy-momentum tensor for the matter Lagrangian Lmat 
is constructed as usual via (I60p . From the Noether theorem we hnd that it 
satisfies the generalized conservation law 


VHij = ^ {4tij - gijt) d"F. 


Here t := g^Hij. The derivation is given in 

To begin with, we recast (j202[l into an equivalent form 


by introducing 


J' = -Ai (S*^' + 

A:=dilogF-^ F}^:=-g^h/4. 


( 202 ) 


(203) 


(204) 


7.3 General equations of motion: arbitrary multipole order 

We now derive the equations of motion for extended test bodies in scalar- 
tensor gravity by making use of the master equations obtained for the general 
case of MAG. In scalar-tensor theory, the hypermomentum is zero = 0. 
Following the general scheme, we introduce the dynamical current which is 
a special case of the MAG current: 




0 

Ykj 


Taking into account the structure of the conservation laws 
we define the material current 

K^=[ 

' -t/4 


(205) 
and (|204p . 

(206) 
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we then recast the system (12031) into the generic conservation law, where we 
now have 


B 


0 


0 

Vjk" 


:k 

°k/ 

0 

0 

-JfT 


In accordance with (j203p and (I204p we now have 


Vjn’^ = Aj5l A = dilogF-\ 


In view of (|205p - (l206p the generalized moments read 


jyi---ynY 

jyi-'-ynYyo 

jyiVi—ynY 


P- 


0 

,yi---yny' 
0 


]^yi---yny yo 






yn 


(207) 

(208) 

(209) 

( 210 ) 
( 211 ) 
( 212 ) 


It is worthwhile to notice that although we formally use a different notation 
for the upper and lower components of the moment (|212ll . they are not 
independent. Recalling (j206l) . we have the obvious relation 

= -^gy,y„^,y^-y^y'y\ ( 213 ) 

which we will take into account when rewriting the equations of motion. 

Now we can make use of the general MAG equations of motion (I119p . 
The first equation is a degenerate version of (|144p since we have = 

0 and qy^---yn-^y<^ybyn = a,re left with the algebraic relation 


j^yi---ynybya _ ^yi---ynyayb 


(214) 


for all n. Hence we find 

( 215 ) 

and the second equation (jl45|) then gives rise to the equations of motion in 
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scalar-tensor theory 

P_pyi-yriya _ _ ^y{yipy2---yn)ya ^j^iyi--yn-i\ya\yn) 

ds^ 


_Ayadyi---yn _ Aya tyi...yn + l 

^ ‘iVn + l^ 

(ky^-y'-+^y'y” + yy” 


-Vy^^y^y-^ky^-y^y'y" - Vyny,y-.y^^^ky^-y^+^y'y” 


+ E 


k=2 


— i — X^^na^y^ / ]^y2---yn)yn+i-yn+kyay' 

V ■'■/ y yn+l-'-Vn + k^ 


^yi-yn+ky'y” _j_ ^y”pyi-y-n+ky'^ 


+(-l)‘7''-.' 


y y yn+i---yn+k 


-K 


,ya 


y y ; 3 /n + l--- 3 /n + fc 


j^yi---yn+ky y — jxy°- 


yn+l---yn+k 


yi •■■yn+k 


^yi--- 


(216) 

The system (I216p is valid up to any multipole order. In the following we 
specialize it to the dipole and the monopole case. 

7.4 Pole-dipole equations of motion in scalar-tensor theory 

For n = 2,1,0 equation (12161) yields 


0 


(217) 

D,b 

ds 

= - Vdc^k^'^^, 

(218) 


= 



- \R\db . 

(219) 


Taking into account that = q, we resolve (|217l) in a standard way to 
find explicitly 

^abc ^ ^a^ch ^ ^c^[ab] ^ ^b^[ac] ^ ^a^[hc] _ (220) 

In view of (I215p . we have in addition 


C = - ^9abk' 


ab 


( 221 ) 

( 222 ) 


Then repeating the same algebra as we did in MAG, we recast the system 
p218l) and p219p into 

1 


^ = ^R\cdV^J'^^ - 

as 2 


DJ' 


ab 


ds 


= _ F~‘^. 


(223) 

(224) 
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Here, following m, we have the total orbital angular moment := 
Whereas the generalized total energy-momentum 4-vector and the general¬ 
ized total angular momentum are introduced by 

;= F" V + (225) 

jab p-iiab_ ^ 226 ) 


7.5 Monopolar equations of motion in scalar-tensor theory 

At the monopolar order we find from eq. (|216l) for n = 1, n = 0: 

- uV, (227) 

- (228) 

Making use of = 0, the first equation yields = 0, hence we have 


0 = 
Dp^ 
ds 




(229) 


with the mass M := v°‘Pa- Substituting (I209D . (|227D and (I229D into 
we find 

D{F-^Mv^] 


ds 

Contracting this with Va, we derive 

D{F-^M) 

ds 

we write (I230p in the final form 




= -iv^V^F-\ 


Dv- - ab,"^bF-^ 


Combining (j222[l with (I227p we find 


^ = - 


V°‘Pa M 


4 4 

Substituting this into (I232p . we obtain 


'ab ah.'^bF 


ds 


= - (g“'’ - vV 


(230) 


(231) 


(232) 


(233) 


(234) 


Quite remarkably, we thus find that the dynamics of an extended test body 
in the monopole approximation is independent of the body’s mass. In case 
of a trivial coupling function F, equation (I234p reproduces the well known 
general relativistic result. 

Interestingly, the mass of a body is not constant. Substituting (|233l) into 
(12311) we can solve the resulting differential equation to find explicitly the 
dependence of mass on the scalar function: M = F^Mq with Mq =const. 
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Table 1: Overview: MAG equations of motion. 

Lagrangian and conservation laws 

L = V{gij,Rijk,Nki^) + F{gij , Rkli^, Nki"-) Lraa.ti'ip'^ 

V,{FA\^) = - 4 * + - Nnk^^^m^) 

V,{F(Sfc^' + = - FA^JiRkim^ - VuNim^) - L^^t^kF 

Equations of motion (any order) 

See equations ()144l) and ()145p . 

Equations of motion (pole-dipole order) 

^ = \R\cdV^J^'^ + Fq<^^<^V‘^Ndch - 

+ 2T(g'='^[“iV'’lcd + - 2^[“V^]E 

j^ah ^ 2p[“^l = -2/i[“''l 

-pa ^ = FiL^’^ + 5’“'') 

Equations of motion (monopole order) 

M^ = - - n“n^) ^ M = p^Va 
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Table 2: Overview: Scalar-tensor equations of motion (Einstein frame). 
Lagrangian and conservation law 

^ ^ + g'hABdiif^djtp^ -2U^ + d'lp, F~‘^gij) 

F = F{ip^) U = U{ip^) -fAB = 1 ab{^^) 

VTjj = {4tij - gijt) d^F 
Equations of motion (any order) 

See equation (I216p . 

Equations of motion (pole-dipole order) 

^ = lR\cdv'’J^'^ - 

- 2^[“V^1 F-* 

Lab ^ 

-pa _ p-ipa _|_ pba^^p-4: jab _ p-ipab 


Equations of motion (monopole order) 
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8 Conclusions 


We have presented a general multipolar framework of covariant test body 
equations of motion for standard metric-affine gravity, as well as its exten¬ 
sions with nonminimal coupling between matter and gravity. Our results 
cover gauge theories of gravity (based on spacetime symmetry groups), and 
various so-called f{R) models (and their generalizations), as well as scalar- 
tensor gravity. 

Our results unify and extend a whole set of works [70l[7Illl2l[53l[73l[7Il 
ESI [Ml Eel EH Ea [78]. In particular they can be viewed as a completion of 
the program initiated in m. in which a noncovariant Papapetrou [2] type 
of approach was used. The general equations of motion (I144D and (I145h 
cover all of the previously reported cases. As demonstrated explicitly, the 
master equation (|119l) allows for a quick adoption to any physical theory, as 
soon as the conservation laws and (multi-)current structure is hxed. Table 
El contains an overview of our main results in the context of metric-affine 
gravity. In particular, the explicit equations of motion at the pole-dipole and 
monopole order are given. We stress that these equations hold for the most 
general case, i.e. including a general nonminimal coupling between matter 
and gravity. The results for standard (minimal) MAG are easily recovered 
by choosing a trivial coupling function F. 

Our analysis reveals how the new geometrical structures in generalized 
theories of gravity couple to matter, which in turn should be used for the 
design of experimental tests of gravity beyond the Einsteinian (purely Rie- 
mannian) geometrical picture. In the case of minimal coupling, we have once 
more confirmed - now in a very general context - that only matter matter 
with microstructure (such as the intrinsic hypermomentum, including spin, 
dilation and shear charges) allows for the detection of post-Riemannian 
structures. However, in gravitational theories with nonminimal coupling, 
there seems to be a loophole which may proof to be interesting for possible 
experiments; i.e. there is an indirect coupling of new geometrical quanti¬ 
ties to regular matter via the nonminimal coupling function F. This may 
be exploited to devise new strategies to detect possible post-Riemannian 
spacetime features in future experiments. 

In addition to the results in MAG, we have explicitly worked out the 
test body equations of motion for a very general class of scalar-tensor grav¬ 
itational theories. Table [2] contains an overview of our main results in the 
context of this theory class. Again the equations of motion at the pole-dipole 
and monopole order for a general coupling function F, which now depends 
on the scalar degrees of freedom, are explicitly given. 

We hope that our covariant unified multipolar framework sheds more 
light on the systematic test of gravitational theories by means of extended 
and microstructured test bodies. We would like to conclude with a statement 
by Einstein m who stressed that 

“[...] the question whether this continuum has a Euclidean, 
Riemannian, or any other structure is a question of physics 
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proper which must be answered by experience, and not a ques¬ 
tion of a convention to be chosen on grounds of mere expediency.” 
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Appendix 


A Conventions &: Symbols 


In the following we summarize our conventions, and collect some frequently 
used formulas. A directory of symbols used throughout the text can be 
found in tables HiaEi 

For an arbitrary /s-tensor f^e symmetrization and antisymmetriza- 

tion are defined by 


T[ai...ak) 


^[ai...afc] 


1 

'kl 


k\ 

E- 

i=i 




1 


k\ 

^^( — 1)1 


I=l 


(235) 

(236) 


where the sum is taken over all possible permutations (symbolically denoted 
by 7r/{ai ... a^}) of its k indices. As is well-known, the number of such per¬ 
mutations is equal to k\. The sign factor depends on whether a permutation 
is even (|7r| = 0) or odd (|7r| = 1). The number of independent compo¬ 
nents of the totally symmetric tensor T'(ai,,,afe) of rank A; in n dimensions is 
equal to the binomial coefficient = (n— 1-|-A:)!/[A:!(n — 1)!], whereas 

the number of independent components of the totally antisymmetric ten¬ 
sor T[ai...as,] of rank A: in n dimensions is equal to the binomial coefficient 
(^) = n!/[A:!(n — A:)!]. For example, for a second rank tensor Tab the sym¬ 
metrization yields a tensor T(^ab) = ^{Tab + Tba) with 10 independent compo¬ 
nents, and the antisymmetrization yields another tensor T^ab] = ^(Tab -Tba) 
with 6 independent components. 

In the derivation of the equations of motion we made use of the bitensor 
formalism, see, e.g., [3 [HI [80] for introductions and references. In particular, 

/ y \ 2 

the world-function is defined as an integral cj(x, y) := i f dr \ over the 
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Table 3: Directory of symbols. 

Symbol 

Explanation 


Geometrical quantities 


s 

Coordinates, proper time 

9ab 

Metric 


Kronecker symbol 

9 

Determinant of the metric 

a 

World-function 


Connection 

Tab^ 

Riemannian connection 

Tab^ 

Transposed connection 

T) d 

^abc 

Curvature 

Qabc 

Nonmetricity 

Qa 

Weyl covector 

n~< c 
^ab 

Torsion 

Kab^ 

Contortion 

Nab^ 

Distortion 

^abt ^ 

Ricci tensor / scalar 


Generators of general coordinate transformations 


Gravitational Lagrangian (density) 


Generalized gravitational field momenta (densities) 


Gravitational hypermomentum (density) 

fr ^ 

Generalized gravitational energy-momentum (den¬ 
sity) 

I 

General action 

L,£ 

Total Lagrangian (density) 


General set of fields 

C“ 

(Generalized) Killing vector 


General dynamical currents (densities) 

J 

9ij 

J J 

Jordan frame metric 

L,£ 

Total Lagrangian (density) in Jordan frame 
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Table 4: Directory of symbols. 

Symbol 

Explanation 


Matter quantities 


General matter field 

-^^mat ; '^mat 

Matter Lagrangian (density) 

sp h cr h 

Canonical energy-momentum (density) of matter 

Aac c 

^ b b 

Canonical hypermomentum (density) of matter 

^ab1 ^ab 

Metrical energy-momentum (density) 

Tab" 

Spin current 

A“ 

Dilation current 


Four velocity 

P 

Pressure (hyperfluid) 

Pa 

Momentum density (hyperfluid) 

Ja^ 

Intrinsic hypermomentum density (hyperfluid) 

F,A 

Coupling function 

pab 

Variation of the coupling function 

ja 

Induced conserved current 


General material currents (densities) 


Integrated moments 

h-,q-,n-,C- 


pab 

Electromagnetic field 


Electric current 

j^ab gab 

Total orbital / spin angular momentum 

pa ^ jab 

Generalized total momentum / angular momentum 

'Y'O'fe 

Total hypermomentum 

M 

Generalized testbody mass 

D 

Intrinsic dilation moment 


Multiplet of scalar fields 

K 

Einstein’s gravitational coupling constant 
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Auxiliary quantities 



Directory of symbols. 

Explanation 

Auxiliary variables (Noether identities) 

Partial derivatives of the total Lagrangian 
Auxiliary Lagrangian density 

Partial derivatives of the coupling function 
General multi-indices 

General functions of external classical fields 
Auxiliary variables (MAG equations of motion) 
Derivative of the coupling function 
Expansion coefficients of the parallel propagator 

Decomposition pieces of the q-moments 
d 

Trace of the Q moment 

General function of scalar fields (in Jordan frame) 

General potential of scalar fields (in Jordan frame) 
Auxiliary variable (Einstein frame) 

Auxiliary variables multipole expansions 

(Partial, Lie) derivative 
Govariant derivative 
Riemannian covariant derivative 
Transposed covariant derivative 
Govariant density derivative 
Riemannian covariant density derivative 

Modified covariant density derivative 
Total derivative 
Variation, substantial variation 
(Bi-)Tensor density 
(Generalized) parallel propagator 

Jordan frame quantity 
Riemannian quantity 
Goincidence limit 
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geodesic curve connecting the spacetime points x and y, where e = ±1 
for timelike/spacelike curves. Note that our curvature conventions differ 
from those in [3 ED]. Indices attached to the world-function always denote 
covariant derivatives, at the given point, i.e. ay := Vycr, hence we do not 
make explicit use of the semicolon in case of the world-function. The parallel 
propagator by g'^x{x, y) allows for the parallel transportation of objects along 
the unique geodesic that links the points x and y. For example, given a 
vector at x, the corresponding vector at y is obtained by means of the 
parallel transport along the geodesic curve as x{x,y)V^. For more 

details see, e.g., 13i or section 5 in |80] . A compact summary of useful 
formulas in the context of the bitensor formalism can also be found in the 
appendices A and B of [26] . 

We start by stating, without proof, the following useful rule for a bitensor 
B with arbitrary indices at different points (here just denoted by dots): 

= [B..,y] + [B..,x] . (237) 

Here a coincidence limit of a bitensor B^^^{x,y) is a tensor 

[B ] = lim B {x,y), (238) 

x^y 

determined at y. Furthermore, we collect the following useful identities: 


^yoyixoy2xi — ^yoyiy2xoxi — 

^xoxiyoyiy2 : 

(239) 

g'^^'^’^ax^ax^ = 2a = g^^^^a 

yi ^y 21 

(240) 

II 

o 

II 

W 

s_, 

II 

a 


(241) 

[' 73 ;ix 2 ] — \^yiy2] ~ 9yiy2: 

['7x11/2] ~ [*71/1x2] = ~9yiy2i 

(242) 

['7x1x23:3] — ['7x1x21/3] ~ ['7x11/21/3] ~ [*^1/11/21/3] ~ 0 , 

(243) 

b"%i]= 5 ^%i, [< 7 "%i;. 2 ] 

II 

0 

to 

II 

0 

(244) 

\g °yr,X 2 X:^ = 2 ^^° y^y^y^- 


(245) 
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B Covariant expansions 


Here we briefly summarize the covariant expansions of the second derivative 
of the world-function, and the derivative of the parallel propagator: 


CJ 


Vo — 

xi — 


a 


yi 


fjVO 

y xi\x2 


„yo 

y Xi-,y2 


G^°X,;X2 

G^^Xr,y. 


9^ X, 


oo ^ 


k=2 


yo — 


oo ^ 

E i!' 9 '" 


k=2 


yiy 2 ---yk+i 


, , , ^2/fc+l I 

•yk+1^ ^ y 


0-y^.. .0-yk+i 


/ 1 ~ oo ^ 

ny' ny” -Ryo , „ _^ 2/0 , „ 

y xiy a;2l 2"^ yy ys^ ^z_^ / yy y 3 -yk +2 

^ k=2 


(246) 

(247) 

) 

(248) 


/ 1 ~ °° T \ 

9 ^^ y’y2y3-yk+2^^^' ■ .(249) 

^ k=2 ^ 

OO ^ 

G^'x,g^''x2 Y. ^^°Y'y"y^...yk+2 ^^^-' ( 250 ) 

fc = l 

OO ^ 

El 5'^V'y2j/3...y.+2'^^"- • (251) 


fc=i 


The coefficients a, (3 ,7 in these expansions are polynomials constructed from 
the Riemann curvature tensor and its covariant derivatives. The first coef¬ 
ficients read as follows: 


r^yo — 

yiy2y3 

(?/2?;3)3/i’ 

(252) 

l^^°yiV2y3 ~ 

(i/2j/3)i/l> 

(253) 

lyyo — 

“ yiy2y3yi ~ 

-V(y2 y3y4.)yii 

(254) 

f^^°yiy2y3yA ~ 

~2^(y2^^°y3y4,)yi > 

(255) 

yyVO — 

I yiy2y3y4 ~ 

-V(y3 R^°\yi\y4)y2- 

(256) 


In addition, we also need the covariant expansion of a vector: 

/ _-|\/c 

^x = X ^yovi-'-Vk (257) 

k=0 


C Explicit form 

Here we make contact with our notation in [5T] to facilitate a direct com¬ 
parison to the results in there. 
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We introduce the auxiliary variables 


^yx-.-Vnyo 

Xq 

:= cj^i- 

• • ay'^gy^ 

q/yi---ynyoy' 

^ XQX' 

:= • 

■ ■ ay^gy° 


(258) 

(259) 


Their derivatives 


^,yi---ynyoy' 


XQX'\Z 




^ o-J/i ... • • • ay'^gy°xo 9 ^'X' 

a=l 


+(T^i ■ ■ ■ a 


Vn 


9^°xo\z9^ 


+ 9^\,9r 


, (260) 


'^(jyi ... (jy<^^ ... ay^gy°^^ 

a=l 

+cr^i ■ ■ ■ gy° xo-^z, 


(261) 


can be straightforwardly evaluated by using the expansions from the previ¬ 
ous appendix. 

In terms of (I258p and (I259p the integrated conservation laws pi27p and 
(11281) take the form: 


I / .,yn+i\uyi---ynyoy' , (c:::xox'x2^y 

^ ^ ^ xox';yn+i^ ^^X2’i 




(262) 


^X2 


XQ 


-R^° 


n///rrfrrf/ 


6 


x'x"x"' 1 r^xQ 


- -Q" 

ir 




+ / ^ / vy-^-rx^yx-. 


X2 


I" 


■ynyo 




(263) 


This form allows for a direct comparison to (29) and (30) in [51] , Explicitly, 
in terms of p258p and p259[l the integrated moments from pi4ip - pi43p are 
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given by: 


pyi---ynyo 


j^y2---yn+iyoyi 


fjy2---yn+iyoyi 


qyi—yn+2yoyiy2 


^y2...yn+iyoyi 


^yi---yn 


(_1)« f (264) 

E(r) 

(-1)^ / (265) 

E(r) 

(_1)^ f q,y 2 -yn+iyoyi (266) 
E(r) 

(_l)n f xl,?/3...2;n+22;o2;i^^^^^2;2^^6xoxiX2^x3^5.^^^(267) 
E(r) 

(_i)« f ^y^-y^+m (208) 

E(r) 

(_l)n I'^yi... (269) 

E(r) 
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